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GENERAL NOMENCLATURE 


I, = dynamic pressure impulse 


D @ displacement 

R # ground range 

x ® scaled ground range 
HOB = height of burst 
s,8 
P 


# scaling factors 


i? 
Po ws, 
a” (ay) wis 


T,(*c) + 273 1/2 aa 4 
Sas 


: C Ga») (y) 
a) 14.7 W 


P,* ambient pressure (lbs/sq in) 


a6 = ambient temperature (degrees C) 

W @ weapon yield (KT) 

&» a> as» b, 2, m, A, B = constants in least squares fits 
x,y - independent and dependent variables in least squares fits 
o ®# standard deviation 

o* ® variance 


N = number of data points 
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r= correlation coefficient 


R = multiple linear correlation coefficient 
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GENERAL NOMENCLATURE (Continued) 


r’, R? s coefficient of determination 


S = sum of squares of deviations of data points from fitted curve 


S' = sum of squares of fractional (relative) deviations of data points from 


fitted curve 


1/2 


Ss ’ 
' s andar 
aren N - number of regression coefficients) ee ee of 


1/2 


Ea Cs) x 100 = root mean square percent error 
SO = side on 


FO = face on 


RO = rear on 


See Appendix, "The Least Squares Procedures", for definitions/defining equa- 
tions of the following specific examples of the above general quantities: 


2 2 2 
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SECTION 1 
(This Section is Unclassified) 


INTRODUCTION 


The objective was to construct improved height of burst contours for 
dynamic pressure impulse. In particular, we wished to extend the set of con- 
tours previously generated (Reference 1) to higher values of dynamic pressure 
impulse and to obtain data points along the contours at additional scaled 
burst heights, i.e., to better determine the contour shapes. To do this we 
analyzed 211 of the data available. This includes dynamic pressure-time wave- 
forms for shots not used in Reference 1 (because no trucks were exposed on these 
shots): TEAPOT Hornet-5, Post-11, and Zucchini-14; PLUMBBOB Franklin-2, Wilson-4, 
Hood-6, Kepler-9, and Owens-10. It also includes data for shots in which tanks 
and howitzers were exposed but no dynamic pressure measurements were made: 
GREENHOUSE Easy-2; TUMBLER-SNAPPER Fox-6 and How-7; UPSHOT-KNOTHOLE Annie-1*, 
Nancy-2, Badger-S, and Simon-7. 


* There was one dynamic pressure measurement on this shot. 
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SECTION 2 
(This Section is Unclassified) 


ANALYSIS OF BLAST MEASUREMENTS 
2.1 Measurements For Eve'.ts Not Previously Considered, 


We consider first the shots not previously considered for which some 
dynamic pressure measurements were available. Some of the waveforms were not 
very good owing tv noise (mixed inextricably with real physical irregularities), 
baseline movement, or gage record cutoff. However, many of the waveforms were 
found useable, being a quality similar to, or not significantly inferior to, 
those used in R:ference 1, General information on all of the nuclear tests of 
concern in this report is given in Table 1. Dynamic pressure impulse, scaled 
and unscaled, is given in Table 2 for each Operation/Event, scaled height of 
burst, and ground range for which we have obtained results. Also given in 
Table 2 is the manner in which we obtained the values of dynamic pressure 
impulse. 


We discuss the data for the shots in the order listed in Section l. 


TEAPOT, Hornet-5. As indicated in Table 1 we used BRL gage data.* 
There were two gage results at each of the three ground ranges. In two instan- 
ces the two results are in good agreement. At the closest-in range, 256 meters, 
the two gages differ by about a factor of two but one of the two appears to be 
correct, the other badly in error. We used the average of the two results for 
the 329 and 460 metre ranges, and the apparently correct results for the 256 
metre ground range. 


TEAPOT, Post-ll. For this shot dynamic pressure waveforms are avail- 
able in Reference 2. We used a planimeter to obtain the area under each curve 
and then used the data reduction procedure described in detail in Reference l, 
A detailed discussion of accuracy of this procedure is also given in Reference 
1. For this shot waveforms are given at four ground tanges in Reference 2. 
The waveforms are not good but are of sufficient quality to be useful. As 
indicated in Table 2 we used the average of BRL's result and our result for 
‘ the two ground ranges at which BRL provided results. Our results and BRL's 
results agreed to within about 3% which is quite good; (the two results are 
based on the same gage data but the data reduction procedures are indeperident; 
we have analyzed errors inherent in data reduction in Reference 1, and, in 
many instances, errors can be quite large). 
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For the two ground ranges at which BRL did not 
reduce the data, we used the result of our own data reduction process. 


a 
a 


TEAPOT, Zucchini-14. For this shot we averaged our result with 
the result of BRL's data reduction process for the three listed ground ranges. 


Me 
bat ad 


* BRL gage data used in this report are published in DNA-TR-85-161. 
Furnished by Mr. J. Keefer, BRL. 
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In each case the difference was between 3 and S%. Results are based on two 
gages at the 610 and 700 metre ground ranges and upon one gage at the 794 
metre range. 


PLUMBBOB, Franklin-2. For this shot we used the BRL result for the 
single ground range for which there appeared to be a good measurement. 


PLUMBBOB, Wilson-4. For this shot BRL also provided values of 
dynamic pressure impulse as listed in Table 2. (We integrated the waveform 
curves as a check, however.) 


PLUMBBOB, Hood-6. Of the 6 ranges at which there were gage results 
the data appeared to be valid at 4 locations. There were two gages at each of 
these locations. We used the BRL results as listed in Table 2 (again after an 
integration check). 


PLUMBBOB, Kepler-9. In this case waveforms obtained at three ground 
ranges were quite poor. In one case there was a large disparity between two 
BRL gage results (almost a factor of 4); in a second case there was almost a 
factor of 2 difference between our gage-reduced result and BRL's result. At 
the third ground range, 762 metres, the waveform was somewhat better. The 
percent difference between our result and BRL's result was rather large, about 
20%; in Table 2 we list the average of the two results. We discarded the 
results at the other two ground ranges. 


PLUMBBOB, Owens-10. For this shot we again used the average of our 
gage-reduced result and the BRL result at two of the three ground ranges, 305 
and 518 metres. At the 305 metre range the results are in reasonably good 
agreement (within 9%); at the 518 metre range the two results differ by a 
factor of 1.4 but the waveforms are quite poor. At the 457 metre range our 
value of dynamic pressure impulse was exactly equal to BRL's value of scaled 
dynamic pressure impulse; the BRL value appeared anomalous on a data plot, 
probably as a result of failure to scale the value obtained. In this case we 
used our gage-reduced result. 


2.2 Analysis of Tank Data. 


* ve ee 


We next consider the set of shots listed in Section 1 on which tanks 
were exposed but there were no dynamic pressure measurements. Our rationale 
here is based on a finding we have discussed in detail in Reference 1: the 
displacement which a vehicle exposed to a blast wave suffers can be used as 
a measure of the dynamic pressure impulse it receives. hat is, the vehicle 
can serve as a gage for dynamic pressure impulse if we are able to "calibrate" 
this "gage". The calibration is the curve of dynamic pressure impulse versus 
displacement. (We are ignoring low yield devices where diffraction effects 
play a role.) We have analyzed the procedure and results for 1/4 ton and 2's 
ton trucks in considerable detail in Reference 1. 
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We would not be able to use the procedure for tanks if we had no dynamic eae 
pressure impulse data for the ground ranges at which tanks were exposed. That hte 
is, we would have no way of obtaining the necessary calibration curve. For- me 
tunately, we have the necessary dynamic pressure impulse data. It is supplied Sh 
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by information from TEAPOT, Met-12 and Apple II-13 and from UPSHOT-KNOTHOLE, 
Annie-1 and Grable-10. The data for the pair of UPSHOT-KNOTHOLE events are 
gage data supplied by BRL; for Grable-10 we also have results based on both 
gage and (truck) displacement from Reference 1, but at larger values of ground 
range than those at which tanks were exposed. The BRL data and the Reference 1 
data are very compatible, i.e., very well represented by a single curve of 
scaled dynamic pressure impulse versus scaled ground range (as shown in Figure 


7, to be discussed later). The results are listed in Table 2. 


For Met-12 we have the results from Reference 1, listed in Table 2, based 
on both gage and (truck) displacement data. (There is a considerably larger 
body of truck displacement data than of tank displacement data so that the 
dynamic pressure impulse - displacement curve, i.e., the calibration curve, is 
better determined for trucks.) 


For Apple II-13 the results listed in Table 2 are taken only partly from 
Reference 1. In this instance, values of dynamic pressure impulse at the 518 
and 625 metre ranges listed in Reference 1 (based only upon gage data) have 
been averaged with BRL gage results, i.e., the two results given equal weight. 
At the 625 metre range this makes a negligible difference while at the 518 
metre range the value of dynamic pressure impulse listed in Reference 1 differs 
from that of Table 2 by about 16%. Finally at the 808 metre range the result 
listed in Table 2 differs from that listed in Reference 1; the latter is in 
error as a check of our (previous) data reduction reveals. 


2.2.1 Dynamic Pressure Impulse Versus Displacement For Tanks. 


We plotted dynamic pressure impulse versus displacement for the 
tanks exposed on events Met-12, Apple II-13, Annie-1 and Grable-10. This 
provides a total of 14 data points* listed in Table 3 along with the tank fae 
displacement data for other shots. Ve 
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There are several variables which cause deviations of the data 
relative to a smooth curve through the points: (1) vehicle orientation with 
respect to the bomb -- 4 data points correspond to side-on orientation, 7 to 
face-on orientation, 1 to rear-on orientation, 1 to face-on 45° orientation, 


and 1 to face-on 3/4 left orientation; (2) there are three different tanks -- we 
M4A3, M24, and M48; (3) there are two surface conditions -- rough sand and may 
fine sand. There obviously are not sufficient data to sort out the effects a) 
of the variables. However, when we plotted the data we found that: (1) Aon 


(initial) orientation appears to make little difference, any effect being 


submerged in effects of the other variables; (we intuitively expect initial 5% 
orientation to make less difference for tanks than for trucks); (2) any eye 
systematic deviation which could be attributed to type of vehicle could as Nos 
well be attributed to inaccuracy in dynamic pressure impulse -~ for example, toes 
the Met-12 points (points B, 6, and 17 in Table 3) are a little high and are wat 
w -q 
* We have not used Point 8 of Table 3, a Smoky-15 point. Placing this point ne 
on Figure 1 shows that there is clearly something wrong with it -- and what oe - 
is wrong with it involves the displacement, not the dynamic pressure impulse, Sears 

even though the latter involves extrapolation on Figure 6. + 
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M48 points while the Apple II-13 points are a little low and pertain to both 
M48's and h24's. (In fact, if we use an average dynamic pressure impulse 
versus scaled ground range curve for the two shots, Met-12 and Apple II-13, 
these deviations disappear.) 


Plots of the data and the fits obtained are shown in Figure 1.* 


We used these data to infer dynamic pressure impulse from the 
measured displacements on events GREENHOUSE Easy-2; UPSHOT-KNOTHOLE Nancy-2, 
Badger-5, and Simon-7; TUMBLER-SNAPPER Fox-6 and How-7. 


We have noted in Reference 1 that, despite rather wide fluctuations 
in displacement owing to the several variable factors mentioned, the central 
curve following the trend of the data points provides a rather good calibration 
which enables us to use a measured displacement to determine the dynamic pres- 
sure impulse to which the vehicle was subjected. Since displacements vary owing 
to uncontrolled factors, the inferred dynamic pressure impulses will exhibit 
considerable dispersion. (However, in many cases described in Reference 1, the 
displacement-inferred values are about as reliable as the gage-inferred values, 
especially when there are two or more values which can be averaged.) 


The calibration curve for tanks is less well determined than that 
for trucks (Reference 1), there being Sewer data points. Also the displacements 
are smaller resulting in greater errors especially from small yield. Nonetheless, 
when the displacement-inferred dynamic pressure impulse values are plotted versus 
scaled ground range (along with a few deta peints for which there also are gage 
data), we see that the curve obtained is reasonably well determined. See Figure 
3. 


We used an eye-drawn curve rather than a proportional fit in Figure lL 
in inferring dynamic pressure impulse from measured displacements for the follow- 
ing reasons: 


(1) point 20, very small displacement, is much less important 
than the other data points in application of this fit, i.e., inference of 
dynamic pressure impuise from measured displacement; 


(2) point 21 has virtually the same displacement as the cen- 
troid so that it can be ignored in drawing a straight line through the data 
points; [the line must pass through the centroid: 


wn D = 2.448, kn Ty = 3.768; 


on th: abscissa and ordinate scales from Figure 1 we see that both fits do 
Paes: Fprough: the porn’ ranch): 11 sSbs! In Ty = 43.295 the lines do not 


pass through the point D = 21.85, T 47.95); 


* Because of the mentioned variables whose effects we cannot disentangle owing 
to the small amount of data, we have numbered the points on Figure 1. Cor- 
responding numbers are listed in Table 3. By comparison one can verify 
that the data do not exhibit any marked effects owing to differences among 
the above variables. 
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(3) the Grable-10 points 21 and 22 are badly inconsistent with 
one another and points 9 and 10 are somewhat inconsistent. Thus our eye-drawn 
curve was a little higher than Fit 1 of Figure 1 at the small displacement end 
and, like Fit 1, passed through the centroid of the data points. The differ- 
ence, however, is nowhere greater than a few percent. We also note from 
Figure 1 that a curve through the data points would do no better than a 
straight line and is therefore not warranted. It is quite possible, however, 
that if the data covered a wider range a straight line would no longer be 


adequate. 
2.3 Discussion Of and Minor Improvements Upon Previously Used Blast Data. 


Most of the dynamic pressure impulse data listed in Table 2 which we 
have not yet discussed is taken from Reference 1. In a few instances a value 
taken from Reference 1 (gage only - no vehicle displacement data available) 
was averaged with a BRL gage-reduced result as shown in the table. In one 
instance, Turk-4, ground range 595 metres, we reduced the gage data even 
though the gage had cut off near the end. We estimated the shape and time of 
pulse completion from other Turk-4 data and believe that the inaccuracy so 
introduced is small — not more than a few percent. Reference to Figure 3 
shows why this is desirable: Turk-4 has the lowest scaled burst height for 
which we have gage data (except for a single Annie-1 point and except for 
surface burst data); its scaled burst height is only moderately greater than 
those of other events shown on Figure 3; the Turk-4 data are consistent with 
the data for the other events shown on Figure 3 and extend to a much greater 
value of scaled ground range. Thus the curve fit of Figure 3 is not diminished 
in reliability by inclusion of Turk-4 data which, however, enable its use over 
a much greater range in constructing height of burst charts, the object of this 


report. 


For events Yuma-4 and Wasp Prime-9 we used our values from Reference 1, 
therein described as "first iteration” results; i.e., when available, truck 
displacement-inferred dynamic pressure impulse values were averaged with gage 
values in obtaining our best estimates of dynamic pressure impulse. (The 
manner of averaging, justification, and discussion of accuracy and reliability 
are considered in detail in Reference 1.) We originally considered Yuma-4 to 
be an ideal event, but we now believe the values achieved by averaging displace- 
ment-inferred dynamic pressure impulse values with the gage values represent an 
improvement over the gage values alone; the results achieved are certainly more 
compatible with the other data in the same scaled height of burst region. See 
Figure 9 (and compare with gage only data listed in Reference 1, page 115). ~ 


For event Encore-9, a near-ideal event, we use the gage results -- even 
though in Reference 1 we sought to improve upon these results by use of 
displacement-inferred dynamic pressure impulse. The reasoning used in Reference 
1 was that since the data scatter on a displacement versus dynamic pressure 
impulse plot was just as great for ideal/near-ideal shots as for non-ideal and 
since the displacement-inferred dynamic pressure impulse for non-ideal shots 
is, on the average, as accurate and reliable as the gage values -- then 
averaging of displacement-inferred with gage-inferred dynamic pressure impulse 
should improve the values for near-ideal events just as it does for non-ideal. 
The points at issue are discussed in detail in Reference 1. One point, however, 
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is that there is a much greater variety of ideal/near-ideal shots than of 
non-ideal, especially regarding the range of weapon yields in the data base. 
When good gage data are available, these data should produce better results 
than can be obtained by averaging in displacement-inferred values -- owing to 
the effects of the several uncontrolled variables previously described. That 
is, while, on the average, displacement-inferred values are reliable, consider- 
able dispersion is to be expected. Compare, for example, Figure 3; here with 
the exception of the Turk-4 data points and a single Annie-1 point, all of the - 
plotted points are based entirely upon displacement-inferred values of dynamic 
pressure impulse. While the central curve following the general trend of the 
data is reasonably well determined, the data point deviations are rather large 
-- considerably larger than the deviations in Figures 4-10. 


Near~ideal shots also have (smail) real differences relative to one 
another. This is shown, for example, by Figure 3.13 of Reference 1 which is 
a plot of scaled dynamic pressure impulse versus scaled ground range. Data 
for the near-ideal (surface burst) shots do not completely coalesce under the 
scaling -- there are clearly small but real differences owing to variable 
factors not controlled in the experiments. In Figures 3-10 the data for the 
events plotted on each figure seem to coalesce quite well with respect to a 
single curve. In Figure 7 some deviation can be seen. In general, deviations 
tend to appear when data for several events are shown on a single plot, when 
there are several data points for each event, and where the data for several 
events overlap, i.e., cover the same domain of the abscissa, scaled ground 
range. The data in Figures 3-10 pertain to a much smaller range of weapon 
yields than do the surface burst data of Reference 1 (Figures 3.13 and 3.15). 
In general, when the data for sever: i shots coalesce we can attach a high 
degree of reliability to the data. 


In constructing height of burst charts in this report we used the results 
described and listed in Table 3. For the surface burst data we used Figure 


3.13 of Reference 1.* 


Finally we analyzed field data for many more shots, but found the data 
unusable. In some instances the waveforms were very bad; in other instances 
the gages were placed at elevations other than 3 feet above the ground. The 
effect of gage elevation is non-negligible. (Different types of gages also 
exhibit somewhat different responses which must be accounted for in order to 
achieve consistent results.) ye 
NOTE: The BRL data used in Table 3 have not been published. 

2.4 Dynamic Pressure Impulse Versus Displacement for Self-Propelled Howitzers. 

Figure 2 is a plot of displacement - dynamic pressure impulse data for 
* For a single (extrapolated) point on the charts for which the scaled dynamic 

pressure impulse is 15 kPa-sec and the scaled height of burst ‘is zero, we 
also used Figure 3.15 of Reference 1 and took account of the difference 
between Figures 3.13 and 3.15 in this region. Figure 3.15 extends to slight- 
ly higher values of dynamic pressure impulse than does Figure 3.13 so that 
less extrapolation is required. 
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self-propelled howitzers. Here we used values of dynamic pressure impulse 
already described. There are only eight data points. These involve 2 vehicles, 
3 orientations, 2 soil types and 7 nuclear events. Three of the data points 

(2 Badger-S and 1 Simon-7) result from shots for which there were no dynamic 
pressure measurements. Therefore, the values we used, taken from Figure 3 and 
from the scaled height of burst chart, Figure 11, are actually based upon the 
displacement-inferred dynamic pressure impulse values for tanks. 


We did not attempt to use the howitzer displacement data to infer values 
of dynamic pressure impulse. To do so, we should proceed in the same manner 
as for the tank data. Referring to Figure 2, we see that of the eight data 
points only five (points 31-35) could be used. No curve of any reasonable 
degree of reliability can be drawn based on these points. 
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SECTION 3 
(This Section is Unclassified) 


| RESULTS AND CONCLUSIONS 


The results described are given in Tables 2 and 3. Least squares fits 

are given in Tables 4 and 5. A description of the methods used in obtaining 
these fits and the definitions and explanations of the several measures of 
goodness of fit listed are given in the Appendix. The simplest measure of 
goodness of fit to understand is the root mcan square percent error (E, and 
Er,). From Table 4 we see that for the displacement - dynamic pressure impulse 

fits, the errors are quite large.* The error in estimating dynamic pres- 
sure impulse from given displacement is fairly large (35-40%) but tolerable 
and as we have seen the values obtained, on the average, are reasonable and of 
considerable value especially when gage data are poor or lacking. If we were, 
however, to attempt to infer displacement from a given dynamic pressure impulse 
the errors would be much larger (150%) -- though again the results would repre- 
sent reasonable approximations to the average displacement to be expected. (The 
reason for the wide disparity between the errors for the two types of inference 
{is seen to be simply due cto the slope of the fitted lines on the log-log plots, 
Figures 1 and 2, or equivalently from the exponents in the fits listed on the 
Figures and in Table 4.) 


| 
| 


The percent errors, Ey_, shown, in Table 5 for the dynamic pressure impulse 
versus scaled ground range fits are quite modest, varying from about 4 to 254, 
These errors are quite comparable to those obtained in Reference 1 for the var- 
ious fits therein. 
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The results described are shown along with the fits, in Figures 1 - 10. 
These results were then used, along with the surface burst data in Reference 1 
(Figure 3.13 and in one instance Figure 3.15) to construct the desired scaled 
dynamic pressure impulse contours as shown in Figures 11 and 12. Also shown 
is the locus of points separating the regular and Mach reflection regions. 


We have drawn the contours to conform as accurately as possible to the 
plotted points while maintaining smoothness and a continuous variation in 
contour shape as we proceed from low to high values of scaled dynamic pressure 
impulse. 


Figure 11 is the better and more reliable of the two charts. The curves 
fit the data points much better than for Figure 12. In the latter case we have 
not attempted to force the curves through the data points as this would lead to 
structure in the contour shapes which is not justified by the data. 


There are various degrees of dust**, the Met-12 data corresponding to the 
heaviest dust case while the Apple II-13 and Bee~6 data represent more moderate 


* The errors here are much larger than in Reference 1 because there are much 
more data for trucks (of concern in Reference 1) than for the tanks and 
howitzers of concern here. 

** We now believe that the degree of "dustiness" is an indicator of the severity 

of the precursor which is related in turn to the temperature of the pre- 
shock thermal layer, 
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dust. This is the principal cause of the difficulty in Figure 12. However, we 
do not believe that separate charts for moderate dust and for heavy dust are 
warranted, This distinction is simply too fine considering the quality of the 


data. 


In conforming the contours to the data points in as reasonable a manner 
as possible, we note that the two charts differ for the higher scaled burst 
heights, i.e., the comparable contours do not coincide, even though the same 
data are used for the two highest scaled burst heights on Figures 11 and 12. 
We cannot regard the contours in this region as well determined; our preference 
in this region is for Figure 11, because of the greater smoothness of the con- 
tours and because we believe the data at the third highest scaled burst height 
are better for Figure ll than for Figure 12. 


Note: In applications of Figures 11 and 12 the scaling factor S, (see 
General Nomenclature) should be used rather than W-/3 for scaling grdund range 
and height of burst; Sj; should be used in conjunction with contours rather than 
w-l/3 for instances in which ambient pressure and temperature are specified and 
differ significantly from the standard values of 14.7 psi (101.4 kPa) and 15°C, 
respectively. 
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SECTION 4 
(This Section is Unclassified) 


RECOMMENDATION 


It is recommended that the Height of Burst Chart, Figure 11, be accepted 
as standard, i.e., as the most accurate obtainable from the totality of 
existent blast wave measurements. Figure 12 which includes Tq for moderate/ 
heavy dust environments is for unique conditions. The description of these 
conditions was discussed in Reference 1. Thus, Figure 12 should be used only 
for similar conditions or be ignored. 
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APPENDIX A 
(This Appendix is Unclassified) 


THE LEAST SQUARES PROCEDURES 


1. INTRODUCTION. 
All of the least squares fits in this report were obtained with the TI-59 
calculator using library programs supplied by Texas Instruments. 
All of the fits were either of the forn 
y=2=mx-e+b (1) 
or y* a) + aX, * anx, (2) 
where y was the natural logarithm of the desired quantity; x was either the 
independent (controlled) variable or its logarithm and the least squares 
variables Xpe X involved only the special case X), "x and xX." x?, 

The least squares equations are obtained by minimizing the sum of the 
squares of the deviations of the data from the curve. If the experimental 
points have a variable scatter for a given small range of x, weights should 
be applied in the procedures so that the least squares equations correspond- 
ing to Equation (1), for example, are obtained by setting the derivatives of 


$= mW [y, - (ax, + »)}* (3) 


with respect to m and b equal to zero, thus obtaining two linear algebraic 
equations which can be solved for m and b. Here S is the weighted sum of the 


squares of the residuals about the fitted curve and 
Wie 07/0," 

where o is a constant (to be determined from the deviations of the data from 

the fitted curve) and a,” is a measure of the expected deviation from the 

true value for an observation Yi (of unit weight). 

In the Texas Instruments programs W, = 1. In individual cases to be 
discussed we will point out that the progrcems used are nonetheless quite 
adequate for our needs. This is mainly due to the fact that, in all cases, 
we minimized the sum of the squares of the logarithm of the desired function. 
This is equivalent to minimizing the sum ot the squares of the percentage 


deviations from the curve since 
a(en y) 2 Y (4) 
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That is, for small deviations from the curve, the fractional error dy is 
equal to A(2n y). For large deviatiuns this is not exact and in our tabula- 


tion of results we show both : 


oe y 
on y, 7 fn Yo ond Haare —) 


where y, 7 the value of y(x;) curresponding to the curve fit. * 
The statistical analog of Equation (4) (Reference 3) is 


Var (fn y) = + Var y (Var = variance) . (4a) 
Y 


Now in all of the cases the values of y covered a large range. If we con- 
sider, for example, a displacement measurement of 100 metres to be in error 
by 2%, the square of its deviation from the true value is 4; if we consider a 
displacement measurement of 1 metre to be in error by 20%, the square of its 
deviation from the true value is 0.04. So even though the percent errors in 
our data are not necessarily uniform (as best these errors are known), assum- 
ing them to be uniform is much closer to reality than any other assumption we 
can make. From the example just given we readily see that the coefficients 
in the least squares fits would be determined almost entirely by the data 
with large values of y it we were to minimize absolute rather than percent 
deviations from the fitted curves. 

2. DYNAMIC PRESSURE IMPULSE VERSUS DISPLACEMENT (a) 

OLSPLACEMENT VERSUS DYNAMIC PRESSURE IMPULSE (b) 

In application of our results we need fits of the data for cases (a) and 
(b), i.e., with each quantity used in the role of independent (controlled) 
and of dependent variable. 

In each subcase of this case the data are fitted very well by straight 
lines on log-log plots. Thus for fitting on the TI-59 calculator the fit is 


of the form y = mx + b (1) a 
with yslnig, x = nD for Case (a) me 
and y=2nDd, x =n Ig for Case (b) ( 
and b2@nsB. (S) ~*q 


* In most cases we omit the subscript i denoting the ith data point; the a 
e 


summations are taken over the data points in all cases so that omission 


of the subscript { will not lead to any confusion. 
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Equation (1) can then also be written [for Case (a) ], I, * Bo" (la) 
with a similar relation for Case (1b). However, after we have fitted Case (a) 
we can calculate (b) directly (without fitting it) provided we recall certain 
data summations from the machine memory. For Equation (1) we obtain from the 


TI-59 fit: 
slope = m = Say xy /N (6) 
=x" + (Ex) °/N 


) eh 2 LY mex 
Yintercept b N (7) 
1/2 
! a,” 
and the correlation coefficient, r =m (=) ‘ (8) 
(oJ 
y 


r* is called the coefficient of determination. The various summations can be 


recalled from the machine memory. a,’ and o are the variances of the 


xX~-array and y-array data and are given by * 


Oe ee ei 


2 1 —\2 
oO, = 7 S(xy - x) (9) a 
2.1 =, 2 
oy =H Ey, - Y) (10) 
with Repu (11) 
and ¥*_q¥ (12) 
N is the number of data points. et, 
. ‘td 
Equations (6) to (8) are not symmetrical in x and y and the algebraic Lay 
inverse of the least squares fit I, = gp" (1a) ne 
iy 4 
l l/a1/m “ae 
is Ds (@) I, but this is not a least 
ie 
squares fit to the data with the roles of Iq and D reversed. Ne 
However, manipulation of Equations (6) to (10) leads to the equation 
Ie PORE TEN, iter Po, as te 
r= [xy - = Exly] ( ee q (ix "]Uzy" - § (ty) 1) (8a) _ 
ak 


cd 


5 
- a 


shen x and y are interchanged. 


ia i 


which is ea in x and y, i.e., unchan- 


ay 
a 


? i “s vire 
A 


* Bar over a quantity indicates an average value. 
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Thus, if we have obtained the least squares fit Ig = BD” we know that 
the correlation ratio r has the same value with the roles of Iq and D 
reversed. It can be. shown (by taking the origin at the data centroid), 
(References 3 and 4), that the least squares fits 


m 
Tq = BD (1a) 
£ 
D2 AI 
q (1b) 
have values of 2 and m which satisfy Emer’, (13) 


Thus, having determined m (and r) in the least squares fit to case la 
we obtain 2 in the least squares fit to Case (1b) from Equation (13). The 
constant A in Equation lb is then determined by the fact that both least 
squares fits (la) and (1b) pass through the centroid of the data distribu- 
tion. (Compare Reference 5.) Thus, from Equation (1b) we have 


inD=Mnarl tl, (1b) 
and in particular Kd = hark &% T (14) 


where from Equations (11) and (12) 
GHD = 4 Ler D) (15) 


tft zie 


and these values are available from the machine memory after running case 
(1a). Thus from Equation (14) we determine fn A, hence A and the least 
square? fit to case (lb) is fully determined. 

Having obtained the least squares fits, the coefficient of determination 
provides a measure of tlie goodness of fit. By manipulation of Equation (3), 
with Wi =z 1, and Equations (6) through (12) it can be shown that 


L-res—y (16) 
No 
Yy 
In this report we are concerned only with positive correlation, so that the 
nearer r is to +1, the smaller is S, the sum of the squares of the deviations 


from the fitted curve, and hence the better the fit to the data. The expres- 


sion for l - r* can also be put in the form 
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2 
Ity - 
Ely =i ¥) 
=, 2 
tty - y) 
thus expressing 1 - r? as a ratio of the sum of the squares of the deviations 


(17) 


2 
l-r # 


| from the curve to the sum of the squares of the deviations from the mean (of 
) the data). The denominator may be regarded as a normalizing factor. It ‘s 
independent of the functional form used for the fitting function. Without 
such a denominator 1 - r2 would tend to increase with the number of data points 
— even if the data were excellent. 
Now letting Sen a and Son D be the sum of the squares of the deviations 


from the fitted curve ies én Ig and fn D, in turn, play the role of dependent 
variable, we have from Equation (16) 


Stn tq , Send ae 


—=— 
Nop, a NOpn D 


(.- r’) = 


"Here Con ae 5 CEen if - _ (en i) 
(19) 
orn D = : (r(2n py? - = (L2n p)7] 
in accord with Equations (9) through (12). Since we have seen that r is 
unchanged when the roles of D and Iq are reversed, Equation (18) shows that 
the sums of the squares of the deviations from the curves are not the same 


for the two corresponding least squares fits (to the same data) but that 


2 
Sen I Sen I 
peal: pepe: | : (20) 
Sen D of 

fn D 


In all cases treated in this report we found S2y, Iq to be considerably 
less than Son D° 

In an application to be made of the results of this case we wish to 
know the standard deviation of the value of the slope of the curve as well as 
the standard deviation of our observation. We use the following notation 


with the subscript c referring to the curve fit in each instance: 
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Orn L..c 2 variance of an observation of fn Ig relative to the 
’ 
q curve; (note that we have assumed the variance of 


the percent error to be independent of the abscissa 


value) 


2 ‘ : 
a, a, ™ variance of slopes of curve fits. 


Then it can be shown that (Reference 6) 


2 /N2* 
Tey be. Sp); ty 


; (21) 
N-2 
Cen D,c : Sen ay 
N-2 
Sen D / 1 - r* o? 
O° = . an D 
£ Noo, I N<-2 oon t 
q q 
(22) 
2 Sen y / Se Ent 
Using Equation(8) twice, first with D, I, being independent and dependent 
variable and then with their roles reversed, we also see that 
: f\ 
he (23) a 
re m * 
m sf 
In the tabular results for the Displacement - Dynamic Pressure Impulse My 
a 


? 
7: 


data curve fits, we list for each subcase the quantities N, Son I? Ten 1c? 
q q’ 


Son D? "Zn D.C? 2, m, Op, 9, in addition to r and the fitting functions. 


To gain further insight into the reliability of the data fits we have 
also tabulated the following additional quantities defined as follows: 


* Here there are N-2 degrees of freedom; the two degrees of freedom lost 
correspond to the number of regression coefficients (two). 
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2 
p-0 
C 
St = £( ) 
D Dc 
(24) 
2 
Lost 
st = I a) 
q qc 
1/2 
| S5 
| %¢* AEP 
(25) 
s: 1/2 


Q 
~ 
a 
ie) 
a 
- 
dl. 
‘ 
—_ 


and finally the root mean square percent error which is 


Ey s (s/n) 3/2 x 100% 
(28) 


Ey = (St my’? x 1008 
q a. 


The o' quantities, Equation (25), are standard deviations of the fractional 
error (percent error apart from a factor.of 100). Thus when the number of 
data points N is fairly lurge (large enough so that the fit is fairly reli- 
able) the values of % ; "test are almost the same as the values of E, and 
EIy respectively, (aside from the factor of 100 expressing the latter quan- 
tities in percent). 

The reason we have chosen to list these additional quantities, in this 
and in other data fits to be discussed shortly, is that in many cases the 
data scatter is rather large (irrespective of the functional form selected 
for fitting the data). This means that a deviation Afn y in fn y for a given 
data point may differ considerably from Ay/y. To illustrate, suppose the 
curve fit and data point values of a displacement are 10 metres and 6 metres. 
The contribution of this point to the sum of the squares of the deviations 
of {nD is then (fn 10- Qn 6)* = 0.261 (irrespective of which value is the 
curve fit value and which is the data point value). 

The contribution to S' is 


D 


2 
19-6)". yy 
(R-) = 9.160 


if the curve fit value is 19 und the data point value is 6; if, however, the 
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curve fit value is 6 and the data point value is 10, the contribution to So 
i ? 

7 AES) "= 0,444. 

Since the least squares curve tends to pass through the data region with 
a fairly uniform distribution of points on either side, on the average the 
curve fit will be below the data point about as often as above it. Averaging 
the above results yields 9.160 - 0. =» 9.302 which is to be compared with 
the value 0.261 for the logarithmic deviation. Had we chosen the values 10 
and 9.8 instead of 10 and 6, the two compared numbers would be virtually 
identical. 

It is inherent in the nature of the data we are dealing with in this 
report that there is large data scatter while the number of data points is 
quite small. We found (Reference 1), using the same least square procedures 
described herein, that, when the scatter of the data is small, Sp i fairly 
close in value to Son D and similarly S; is fair'y close to Son I.* Usually 
Sen p< S5 and Seon I < ata The situation is sometimes reversed, ‘however, 
usually because one or two data points with large deviations from the curve 
fit are below the curve. Conversely, for a case in which one or two data 
points are far above the curvé it can happen that sy >> Son y since there is 
no limit to the contribution such a point can make to sy y here indicating 
any quantity whose 2n has been fitted. (For a data point below the curve fit 


the maxisum contribution to sy is 
2 2 


The results shown in Table 4 are consistent with these statements taken from 


Refererce 1. [The large effect which a single data point can have may be 
seen by calculating D-D., é 
Cae 

c 


for point 31 of Figure 2. The value, 18.52, is 


almost as large as S) (= 19.06) for the 8 data points of Figure 2.] 

Finally, quantities such as oF Cc and Er_ are directly associated with 
the plotted data and are easily visualized. Since they pertain to fractional 
or percent errors they apply equally well at all points along the curve fit 
(although the curve fit has greater predictive accuracy in the vicinity of 
the data centroid than toward the ends of the data range). On the other 
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hand Ten I refers to the natural logarithm of the plotted data (ordinate) 
q 


and the latter may be quite small or large (or even negative), so that while 
Stn Ig is a standard measure of the data scatter with respect to the fitted 
curve, it is not as simple to interpret as are the additional quantities 
which we have tabulated. 
In sum, the tabulated measures of goodness of fit provide some insight 

into the assessment of the reliabilicy and utility of the data and the 
curve fits and help with the interpretation of results. 

The quantities described here are also tabulated for various curve fits 
of data other than displacement-dynamic pressure impulse data. We discuss 
the least squares treatment of the remaining data in the following section. 


3. ADDITIONAL LEAST SQUARES FITS. 
Scaled Dynamic Pressure Impulse versus Scaled Ground Range: 
The fitted curve is given by 


| 
| 
| 


y= a+ a,x + a5x” (27) 
where x = scaled ground range, y = én S1,- 
In running this case on the TI-59 calculator we feed in the natural 

logarithm of each data point. Again, the machine generates a fit and all of 
the relevant sumaistions over the data are available from the machine memory. 
For this type of fit (trivariate), however, the machine provides a quantity 

R2 rather than r as a measure of the goodness of fit. R is the multiple 
linear correlation coefficient between y and the other least squares variables 


xX, =x and xX. 5 x; R? is called the coefficient of determination (see Ref- 
erence 6) and is given by 
Toy fg Se 28 o5 Fat 
a? » i y2 12 “yl y2 (28) 
1 - r? 
12 
The subscripts 1 and 2 refer to x) and Xi Tote TY2? and T12 are given by 
equations analogous to Equation (8a). 
‘ 2 2 -1/2 
4. 
ry, * [Exy - ExEY (ee ; GA itty? - Gp"y) (29) 
27 ey: 
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2 2,2 2. 
ryg = Cixty - BEE) (tex4 - GeV iery? . Ga"s) 
(29) 


“1/2 
2 2 9 
yp = tex? - Ey (tx? - Ga strt . GO 3) 


N 


that the least squares variables x) and X, as used here are x and x’, respec- 
tively.) The quantities r 1 Tye" and 15 are the coefficients of correlation 
(also called simple correlation or zero-order coefficients) between y and Xp 
y and x,, and % and x», respectively. The important point here is that 


(The superficial lack of symmetry between TH and To results from the fact 


> Elyy,)? 
(1 - R°) = — (17a) 
E(y-y) 
so that R? provides a measure of che goodness of fit precisely similar to that 


provided by r? in the previous discussion: R? is a measure of the closeness 


of fit of the regression plane (in fn y, Xp» Xy space)* to the data points. 
There are two further changes: replacement of N-2 by N-3 in Equations 
(21) and (25) and replacement of I. by S,I 


q q 


2 
Q 2S N-3 
tn S,1,,C ~ "en si1,/ 


(21a) 
of = N-3 
én D,C an D 
and 
1/2 
1c s (Si/N-3) 
25a) 
1/2 ( 
fof a (S2 /N-3) 
$,140C Sit, 
Referring to the results listed in Table 5 we see that usually 
Sen s.t <g: (and in fact this relation holds for all of the cases where 
ivq S,I 
i'q 


the number of data points > 11). The two values are usually quite close, 
however, as might be sxpected, since the data dispersion is fairly modest. 


* Xx; * % = scaled ground range KX, =X 
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